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ABSTRACT 

The  effect  of  small  amplitude,  time-periodic,  freestream  disturbances  on  an 
otherwise  steady  axisymmetric  boundary  layer  on  a  circular  cylinder  is  considered. 
Numerical  solutions  of  the  problem  are  presented,  and  an  asymptotic  solution, 
valid  far  downstream  along  the  axis  of  the  cylinder  is  detailed.  Particular 
emphasis  is  placed  on  the  unsteady  eigensolutions  that  occur  far  downstream,  which 
turn  out  to  be  very  different  from  the  analogous  planar  eigensolutions.  These 
axisymmetric  eigensolutions  are  computed  numerically  and  also  are  described  by 
asymptotic  analyses  valid  for  low  and  high  frequencies  of  oscillation. 
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1,.  Introduction 

The  effect  of  time-periodic  disturbances  in  the  frccstrcam  of  an 
otherwise  steady  boundary  layer  has  received  considerable 
attention  over  ( he  years.  This  work  was  initiated  by  Lighthill  (1954), 
who  considered  t  tic  flow  past  a  semi  -  inf  ini tc  flat  plate,  with  a  small 
amplitude,  time-periodic  frccstrcam  disturbance,  and  obtained  solutions 
close  to  and  far  from  the  leading  edge.  This  work  was  later  extended 
by  Roll  arid  Rosenweig  (1960),  Lain  and  Roll  (1960)  and  Ackcrberg  and 
Phillips  (1972).  Of  particular  interest  arc  the  unsteady  cigcnsolut ions 
that  form  part  of  the  far-downstream  flow.  One  set  of  these  was 
studied  by  Lam  and  Roll  (1960),  Ackcrberg  and  Phillips  (1972)  and 
Goldstein  (1983)  and  has  exponentially  decaying  solution  downstream  (see 
(8.1)  below),  with  the  feature  of  dccrcas ing  decay  rate  with  increasing 
order;  these  cigcnsolut ions  arc  determined  primarily  by  conditions  close 
to  the  wall.  A  second  set  of  cigcnsolut ions  was  constructed  by  Brown 
and  Stcwarlson  ( 1973a, b)  and  has  the  feature  or  increasing  decay  rale  with 
increasing  order:  these  cigcnsolut ions  arc  determined  from  conditions  far 
away  from  the  wall,  in  the  outer  reaches  oT  the  boundary  layer. 

Indeed,  these  seemingly  diverse  characteristics  of  the  cigcnsolut ions 
have  been  the  subject  of  some  controversy  over  the  years.  However, 

Goldstein  cl  ai  (1983)  include  a  quite  detailed  discussion  of  this 
dichotomy;  briefly,  these  authors  expound  the  argument  that  the  two  sets 
of  cigcnsolut ions  arc  in  fact,  equivalent,  but  with  the  Brown  and  Stcwarlson 
(1973  a,b)  expansions  being  valid  at  much  longer  distances 
(0(ln  j  x)i  »  1)  downstream,  than  the  Lam  and  Roll  (I960)  cigcnsolut ions 
(which  arc  valid  for  0(x)  »  1).  Further,  Goldstein  ct  al  (1983)  point 
out  that  as  the  order  of  the  Lam  and  Roll  (I960)  cigcnsolut ions  increases, 
the  asymptotic  behaviour  of  the  (inner)  solution  is  likely  to  be  achieved 
al  progressively  larger  values  of  x,  since,  for  x  »  1,  the  scale  of  the 
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region  associated  with  the  cigensolutions  moves  away  from  the  wall  with 
increasing  order.  This,  in  some  ways  is  not  inconsistent  with  the  fact 
that  the  Brown  and  Stcwartson  (1973a, b)  eigensolutions  are  centered  at  the 
outer  edge  of  the  steady  boundary  layer.  Goldstein  et  al  (1983)  also 
conclude,  using  these  arguments,  that  the  limit  as  x  — *  «  and  the  limit 
as  n  — *  «  (where  n  is  the  order  of  the  eigensolut ion)  cannot  be 
interchanged.  However,  and  significantly,  Goldstein  (1983)  went  on  to 
illustrate  the  physical  importance  of  the  Lam  and  Rott  (1960) 
eigensolutions,  by  showing  how  these  develop,  far  downstream,  into  unstable 
Tol lmien-Schl icht ing  waves. 

The  problem  of  ''order-one"  unsteady,  frccstrcam  disturbances  (but 
such  that  the  freest  ream  docs  not  reverse  direction)  has  been  Considered 
by  a  number  of  authors.  Pcdlcy  (1972)  considered  this  problem, 
asymptotically  close  to  and  far  from  the  leading  edge,  whilst  Phillips 
and  Ackerberg  (1973)  presented  numerical  solutions  to  the  problem  for 
locations  from  the  leading  edge  to  far  downstream,  their  method  being 
based  on  a  time-marching  scheme.  More  recently,  Duck  (1989)  presented 
a  new  numerical  method  to  tackle  this  problem,  based  on  a  spectral 
treatment  in  time,  and  a  spatial  finite-difference  scheme,  which  properly 
takes  into  account  regions  of  reversed  flow  that  inevitably  occur. 

The  problem  of  steady  flow  along  a  circular  cylinder  (in  particular 
far  downstream  along  the  axis  of  the  cylinder)  is  itself  interesting,  partly 
because  it  is  so  very  different  in  nature  from  that  of  planar  (i.e.  Blasius 
type)  flow.  Early  investigations  of  this  problem  include  the  work  of 
Glauert  and  Lighthiil  (1955)  and  Stcwartson  (1955),  whilst  Bush  (1976)  has 
presented  a  more  modern  approach  to  the  problem.  Notably,  in  the  far 
downstream  limit,  the  problem  becomes  double  structured,  with  an  inner 
layer  (comparable  in  thickness  with  the  radius  of  the  body)  which  is 
predominantly  viscous  in  nature,  and  an  outer  layer  (much  larger  than  the 
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radius  of  the  body)  which  is  a  region  of  predominantly  uniform  flow  (see 
Section  4  for  fuller  details). 

In  this  paper  we  investigate  the  effect  of  small  amplitude,  time- 
periodic,  frccstrcam  disturbances  on  the  axisymmctric  boundary  layer  on 
a  circular  cylinder.  Particular  emphasis  is  placed  on  the  cigcnsolut ions 
relevant  to  (he  fa r -downs t ream  flow,  which  turn  out  to  be  markedly 
different  from  the  analogous  planar  cigcnsolut ions  of  Lam  and  Rott  (1960), 
and  possess  some  interesting  properties.  Further,  since  an  additional 
lengthscalc  is  present  in  the  problem  (i.c.  the  body  radius),  a  second 
non-dimensional  parameter  (in  addition  to  the  Reynolds  number)  is  present, 
and  we  arc  able  to  exploit  this  parameter  from  an  asymptotic  point  of  view. 

The  layout  of  the  paper  is  as  follows.  In  Section  2  the  problem  is 
formulated,  and  in  Section  3  a  fully  numerical  finite-difference  scheme 
for  the  steady  and  unsteady  problem  is  described,  and  results  for  the  wall 
shears  arc  presented,  for  axial  locations  from  the  leading  edge  to  far 
downstream.  In  Section  4  the  development  of  the  (inhomogeneous)  component 
of  the  flow  is  described.  In  Section  5  the  presence  of  cigcnsolut ions  far 
downstream  is  elucidated,  the  cigcnproblcm  is  formulated  and  expanded  in  the 
form  of  an  asymptotic  series.  In  Section  6  numerical  solutions  of  the 
(leading-order)  cigcnproblcm  arc  described,  whilst  in  Section  7  the 
cigcnproblcm  is  considered  in  the  asymptotic  limits  of  high  and  low 
frccstrcam  oscillation.  In  Section  8  the  conclusions  of  the  paper  arc 
presented . 
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2.  Formulation 


Wc  introduce  a  cylindrical  polar  coordinate  system  (ar,  9,  az),  where 
a  is  the  radius  of  the  body  (assumed  constant),  and  the  z  axis  lies 
along  the  axis  of  the  body,  with  z  =  0  corresponding  to  the  tip  of 
the  body. 

Suppose  that  the  fluid  is  incompressible  and  of  kinematic  viscosity 
v,  the  frccstrcam  velocity  is  taken  to  be  purely  in  the  z  direction, 
and  of  the  form  W„(l+8  cos  oil*),  where  W„,  5  and  0)  arc  constants, 
with  8  «  1.  Note  that  although  in  all  the  ensuing  analysis  wc  shall 
confine  our  attention  exclusively  to  frccstrcam  velocities  of  the  above 
form,  it  is  relatively  straightforward  to  extend  our  ideas  to  other 
spatial  and  temporal  (periodic)  variations. 

The  velocity  field  is  written  as  (u.O.w),  and  non-dimensional  lime 
as  t  =  mt*.  Further,  it  is  assumed  that  u,w,  and  indeed  the  entire 
solution  is  independent  of  0,  implying  axial  symmetry. 

In  this  problem  there  arc  two  fundamental  non-dimensional  parameters, 
namely  a  Reynolds  number  based  on  cylinder  radius 


which  will  be  assumed  to  be  large  throughout  this  paper,  together  with 
a  frequency  parameter 

P  =  ~  •  (2.2) 

Wit  2 

The  usage  of  the  boundary- layer  approximation  requires  (hat 

Z  =  R-'z  (2.3) 

is  the  key  axial  lenglhscalc,  and 

U  =  Ru  (2.4) 

is  the  important  order -one  radial  velocity  scale.  The  boundary -layer 
equations  then  become  (to  leading  order) 
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J.  +  +  =  ^  +  +  (2.5) 

P  dl  dZ  dr  dr 2  r  dr  p  dt 

together  wi th 

1Z  (rw)  +  I?  (rU)  =  ()-  (2*6) 

Since  it  will  be  assumed  that  8  «  1,  the  unsteady  component  of  the  flow 
may  be  taken  to  be  a  small  perturbation  about  the  steady  solution 
(a  similar  treatment  has  been  used  in  many  of  the  related  planar  studies 
cited  in  the  previous  section,  for  example  Lam  and  Ron  1960,  Lighthill 
1954,  Ackcrbcrg  and  Phillips  1972).  Specifically 


U(r 

,Z,l) 

=  U0(r,z) 

+  8 

Re 

{U(r,Z)e'M 

+  0(82), 

(2.7) 

w(r 

\Z,t) 

=  w0(r,Z) 

+  8 

Re 

{w(r,Z)c' 

+  0(82) • 

(2.8) 

steady  component  of 

the 

solution  is 

described  by 

WQ 

+ 

d2wp 

+  . 

1  dwQ 

(2.9) 

dZ 

dr 

dr 2 

r  dr 

d 

dl 

(r  w0) 

+  Jr  = 

0. 

(2.10) 

W()< 

:  r=l )  = 

=  U0(r=1 ) 

=  0, 

w0 

— »  1 

as  r  — * 

00 , 

(2.11) 

whilst  the  unsteady  perturbation  to  this  flow  is  given  by 

iw  ,  dw  ,  ~  dwn  ,  dw  ,  n  dwn  d2w  1  dw  i 

—  +  w0  —  +  w  — U  +  Uq  —  +  U  — “  =  -  + - +  ~, 

P  dZ  dZ  dr  dr  dr 2  r  dr  p 

(2.12) 

32  (rw)  +  ^  (rO)  =  0,  (2.13) 

subject  to 

w( r=l )  =  U(r=l)  =  0, 

w  — »  1  as  r  — »  ».  (2.14) 

To  close  the  problem  we  further  suppose  that  as  Z  — »  0,  planar 
conditions  prevail,  with  the  boundary- 1  aye r  thickness  becoming 
negligible  in  this  limit.  A  similar  procedure  was  followed  by  Seban  and 
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Bond  (1951)  and  was  further  utilised  in  a  related  problem  by  Duck  and 
Bodonyi  (1986).  The  (steady)  system  (2.9)  -  (2.11)  then  reduces  to  the 
Blasius  (planar)  problem  as  Z  — ♦  0.  with  corrections  due  to  curvature 
effects  given  by  Seban  and  Bond  (1951).  As  Z  — ♦  «•,  the  far  downstream, 
double-structured  solution  of  Glaucrt  and  Lighthill  (1955), 

Stewartson  (1955)  and  Bush  (1976)  emerges  from  this  system. 

Regarding  the  unsteady  system  (2.12)  -  (2.14),  this  becomes  quasi-steady 
in  form  as  Z  — ♦  0,  with  the  time  derivative  term  vanishing  in  this  limit. 

In  the  following  section  fully  numerical  solutions  to  both  the  steady 
and  unsteady  system  and  considered,  and  in  the  later  sections  of  this  paper 
the  far-downstream  behaviour  of  the  unsteady  component  of  the  flow  is 
investigated  in  some  detail. 
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In  this  section  we  consider  fully  numerical  solutions  to  systems 
(2.9)-(2. 1 1 )  and  (2. 12)-(2. 14) . 

Two  strcamfunciions  arc  introduced,  one  for  the  steady  component  of 
the  flow,  the  other  for  the  unsteady  component,  vi/.  4*  and  \jr,  respectively 


given  by 

u<>  -7  S'  D  *  7  S' 

IS- 

The  problem  determining  4*  and  y  is  then 
r  rz 


(3.1) 


(3.2) 


-  7  fr  »rZ  '  fz[frr  '  72  fr]  • 


(3.3) 


»rrr  '  ^  -  -  [»r  VrZ  +  Vr  frz] 


'  [  ?"■  H  ]  ’  «  [  V  ’  %] 


.  |  .  JI 

+  '  t  V 

with  *P  =  *Fr  =  tjr  =  yr  =  0  on  r  =  1 , 
¥r,  5j/r  — »  r  as  r  —*  <*>. 


(3.4) 


(3.5) 

(3.6) 


Anticipating  a  Blasius- type  solution  as  Z  — >  0,  the  problem  for 


0  <  Z  <  1  was  cast  in  terms  of 
r\  =  (r-l)/Zi,  C  =  Zi, 


(3.7) 


as  the  independent  variables,  with  the  dependent  variables  taken  as 


Fo  and  F,  where 


*  =  C  Fo  (i.C).  V  =  CF(n.C) 


(3.8) 


For  Z>  1,  T(r1Z)  and  y( r ,Z)  were  treated  as  the  unknown  variables. 
In  both  Z  <,  1  and  Z  >  1  the  systems  were  written  as  a  system  of 
first  order  equations  in  r  (or  q).  Having  solved  the  problem  for 
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Z  =  0,  a  Crank-Nicol son  procedure  In  Z  (or  £)  was  adopted.  Overall, 
the  numerical  di ffcrcncing  scheme  was  based  on  that  of  Keller  and  Ccbcci 
(1971).  At  each  Z  (or  £)  station  first  the  steady  system  was 
computed,  with  Newton  iteration  being  used  to  treat  the  non-linearity 
in  the  problem.  Once  convergence  was  achieved,  the  (linear)  unsteady 
system  was  then  computed  in  a  straightforward  manner. 

Results  for  Yrr  I  r_j  (essentially  the  steady  component  of  wall  shear) 
along  the  cylinder  arc  shown  in  Pig.l.  This  illustrates  the  (Blasius-lypc) 
singularity  as  Z  — *  0,  together  with  a  monotonic  decline  as  Z  increases. 

Figure  2  shows  the  results  for  the  real  and  imaginary  components  of 
¥rr  lf=i  (essentially  the  unsteady  component  of  wall  shear)  for  p  =  0.25. 
This  shows  how  the  real  component  exhibits  an  inverse  square  singularity 
as  Z  — »  0  (in  line  with  that  of  trr  I  . )  whilst  the  imaginary 
component  drops  to  zero  at  the  leading  edge.  This  occurs  because  as 
stated  previously,  as  Z  — •  0,  the  system  determining  F(q,£)  becomes 
quasi -steady ,  with  the  unsteady  velocity  perturbation  moving  entirely  in 
phase  across  the  boundary  layer.  For  Z£  1,  both  the  real  and  imaginary 
components  rapidly  approach  constant  values.  This  aspect  is  dealt  with 
in  the  following  section. 

Figures  3.4  and  5  show  the  corresponding  distributions  for 
p  =  1,2  and  5  respectively,  all  of  which  exhibit  similar  qualitative 
features  to  the  p  =  0.25  results,  although  the  asymptotic  amplitude 
of  Vrrlr_ i  is  seen  to  diminish  as  p  increases.  In  the  following 
section  the  asymptotic  form  of  the  flow  structure,  far  downstream  of 
the  leading  edge  is  considered. 
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In  this  section  we  investigate  the  2  »  1  solution  for  the 
(unsteady)  system  (2. 1 2) - (2 . 14) .  It  was  shown  by  Glaucrl  and  Lighthill 
(1955),  Stewartson  (1955)  and  Bush  (1976)  that  the  steady  solution 
obtained  from  (2.9)-(2.ll)  divides  into  two  layers  far  downstream. 
Specifically,  for  r  =0(1)  it  was  shown  that 

V  =  l  en+1  Y  (r) 
n=() 

«  n+2 

+  *£  t|n(r) 
n =0 


whe  re  e  = 


l  og  Z’ 


(4.1) 


(4.2) 


and  where 


*0n<D  =  *0n  Ur2  log  r  -  ^  r2  +  ^ 
K00  =  1.  K01  =  (|  -  log  2). 

7  =  0.5772...,  (Eulers  constant), 
■ho  =  K10  li  r2  log  r  •  j.r2  +  J|, 

K10  =  l  • 


(4.3) 


(4.4) 


(Note  that  in  Stewartson  1955,  the  last  term  in  his  equation  (3.20)  should 
be  a  logarithmically  squared  term,  and  not  as  shown).  It  is  also  found 

r  r 

Ti«r>  '  -K00'  r  f  {  jp  I  Yk  fao  Jr  (  Y1  j  ]dl> 


+  K()()r  log  r  j  ^4-Jp  dr  -*■  K, ,  r  log  r  ,  (4.5) 


J  Uf  ** 
1  T00 


implying  that  for  r  »  1 , 
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(4.6) 


oo 

~  l  l  ajn  rj  (log  r)2*n  . 
j=0  n=0 

Consider  now  the  outer  layer,  wherein 

n  =  r/Zi  =  0(1)  (4.7) 

(consistent  with  (3.3)),  wherein 

V  =  Z  {  %,0i)  +  e  Void!)  +  »(e2)  } 

+  e2  f10(n)  +  0(e3).  (4.8) 

A 

it  is  found  f()()(ti)  =  J  q2, 

^oidi)  =  i  ti2  J  ^  dn 

oo 

+  cln2  ■  I. 

»IOOI)  -  H'oiCn)  •  t-  (4-9) 

For  comparison  with  our  fully  numerical  results,  asymptotic 
approximations  to  the  basic  flow  determined  Troni 
Yrr  lr=,  *  «  foorr  («<)  +  e2  foirr  (r*D 

=  e  +  K01  e2  (4.10) 

are  shown  on  Fig.1  as  a  broken  line. 

Now  consider  the  Z  >  1  solution  to  the  system  given  by 
(2. 12)-(2. 14) .  This  turns  out  to  be  quite  straightforward.  Consider 
first  (and  most  importantly)  the  radial  scale  r  =  0(1);  then  due  to  the 
smallness  of  r. ,  w  is  expected  to  develop  as 
w(r,Z)  =  wq( r)  +  0(c), 
where  wq  is  to  be  determined  from 

d2*0  |  1  ^0  j»0-  .  lt 

9r2  r  dr  p  P 

the  appropriate  solution  of  which  is  simply 


(4.11) 

(4.12) 
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wo(r)  =  1 


"0<2)[/rl 


(4.13) 


i.c.  the  axi symmet ric  Stokes  shear-wave  solution,  where  H()(2)(z)  denotes 
the  second  Hankcl  function,  of  order  zero  and  argument  z.  Note  also  that 
as  P  — »  0,  the  planar  Stokes  shear  solution  is  retrieved  (in  accord  with 
the  work  of  Ackcrbcrg  and  Phillips  1972).  In  this  limit,  a  thin  Stokes 
layer  forms  on  the  surface  of  the  body  and  consequently  curvature 
effects  become  less  important. 

It  is  a  routine  matter  to  continue  this  solution  to  higher  orders 
of  e;  however  little  additional  insight  is  gleaned  from  this,  and  instead 
we  go  on  to  consider  (briefly)  the  outer  layer,  where  q  =  0(1)  (see 


(4.3)). 


Writing 

A  A 

w  =  w()(q)  +  e  wj (q)  +  ()(e2),  (4.14) 

then 

w0(tl)  =  1.  (4.15) 

W1 (q)  =  w2(q)  =  . .  =  0.  (4.16) 

In  fact  the  correction  to  wq(q)  can  only  be  algebraically  small  in 


Results  obtained  using  this  asymptotic  structure  (in  particular 
(4.8))  arc  shown  for  comparison  with  the  fully  numerical  results  as 
broken  lines  on  Pigs.  2-5;  the  agreement  is  seen  to  be  satisfactory. 

However,  since  the  Z»  I  structure  detailed  above  is  obtained 
without  any  recourse  to  upstream  conditions,  there  must  be  a  further 
element  to  the  downstream  flow,  not  reflected  in  the  above  analysis 
(see  also  the  comments  of  Ackcrbcrg  and  Phillips  1972).  This 
arises  from  eigenfunctions  of  the  system  (2 . 12)- (2. 14) .  This  aspect 
is  investigated  next,  in  some  detail. 
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Here  the  form  of  (exponentially  small)  cigensolut ions 
as  Z  — »  «•  is  sought.  Specifically,  we  investigate  cigensolut ions  of 

(3.4) ,  with  the  basic  flow  described  by  Section  4. 

As  a  first  approximation  to  the  form  of  these  cigensolut ions , 
consider  the  scale  r  =  0(1),  and  suppose  ¥  in  (3.4)  is  replaced  by 
e¥00(r),  and  terms  O(^)  and  smaller  arc  neglected.  This  yields 

Vrrr  -  ^  -  f  %)r  »rZ 

.  e  lfz  |  Umi  .  iBjl  |  .  o.  (5.1) 

Assuming  a  solution  for  \j?  by  separation  of  variables,  namely 

V=l'(Z)v(r),  (5.2) 

then 

Vrrr-^  +  ?f-  if 

*  e  T2  {  T  f00r  Vr  +  V  I  ^  U 

=0.  (5.3) 

A  solution  of  the  assumed  form  is  possible  only  if 

rz  +  £f  =  o,  (5.4) 

where  A  is  a  constant.  Recalling  the  definition  of  e  in  (4.2), 

(5.4)  integrates  to  give 

f(Z)  =  cxp{-£  |Z  log  Z-Z|} 


Here  it  is  required  that  Rc(A)  >0  to  ensure  decay  as  Z  — ♦  »,  and 
the  arbitrary  multiplicative  constant  in  y(r)  has  been  included. 

However  (5.2)  and  (5.5)  arc  correct  only  to  leading  order  in  e 
and  Z.  It  turns  out  the  form  of  y  required  for  r  =  0(1) 
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V  =  h(Z)  f(Z)  Zp  (log  Z)q  {Voo(r) 
+  e  Yoi(r)  +  ()(e2) 

+  7  le  Vio(r)  +  ()(r.2)  | 

+  0( I /Z2) } , 


(5.6) 


where  f(Z)  is  given  by  (5.5)  and  h(Z)  is  smaller  than  any  power  of 
log  Z.  Further  it  is  found  necessary  to  expand  A  itself  in  terms  of 


ascending  powers  of  e,  viz 

A  =  Aq  +  eAj  +  e2  A2  +  0(e3). 


(5.7) 


p  and  q  arc  constants  to  be  determined  at  some  later  stage.  In  view 
of  our  comments  regarding  Rc(A),  then  Rc(Aq)  >  0.  The  form  of  (5.6) 
and  (5.7)  is  necessitated  because  of  the  series  development  of  the 
basic  flow  in  powers  of  e  and  l/Z,  and  is  found  to  be  essential  for 
solubility  at  higher  orders  of  the  solution.  (Indeed  Goldstein  1983 
pointed  out  the  omission  of  algebraic  terms  in  the  strcamwisc  development 
of  the  planar  cigcnsolul ions  in  the  work  of  Ackcrbcrg  and  Phillips  1972, 
which  contained  only  the  exponential  development  of  the  flow). 

Substitution  of  (5.6)  and  the  results  of  Section  4  into  (3.4),  and 
.  M  AnZ 

taking  terms  o|h(Z)  Z  c  Zp  (log  Z)(*j  yields  the  following 

equation  for  Vqq 


U%>)  =  0. 


(5.8) 


whe  re 


L(*00)  *  Yo<>'  "  -  +  WO'  |  y2  ■  p  *  r°  Y00] 

-*ow»)iy-  y. 

(5.9) 

Recalling  the  form  of  T(j(),  given  in  (4.3),  then 

Mvoo)  =  Voo  ”  ■  ^  +  ¥00 ’  [  72'  P  +  A()  l0«  r] 

-  A()  soil  =  0.  (5.10) 
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The  boundary  conditions  to  be  applied  to  this  system  arc  those  of 


no-slip  and  impermeability  on  r  =  1,  i.c. 
YOO  (r=D  =  Yoo'  (r=0  =  0. 


(5.11) 


whilst  as  r  — >  »  YOO  should  noi  be  exponentially  large.  To  be 

more  precise  on  this  last  point,  the  three  linearly  independent  solutions  to 


(5.10)  in  this  limit  take  the  form 
V()0A  -  A()()A  { 1  og  r  -  jtt~  } , 


(5-12) 


ij  I A()  log  r  |1  dr 
-  ^3/4  •  1 


(5.13) 


- i |  (A()  log  r  ll  dr 


log  r)3/4 


(5.14) 


Clearly  either  one  of  (5.13)  or  (5.14)  is  inadmissible  (if  Aq  is  complex) 


due  to  the  r  »  1  condition,  and  so 

V00  =  A()()  [log  r  -  i-  +  —1 - 

PA0  rzA()log  r 

+  0(  — - - J  ] 

r2(log  r)^  J 


(5.15) 


in  this  limit,  where  Aqo  is  an  arbitary  constant  (amplitude).  The 
system  (5.10)  ,  (5.11)  and  (5.15)  represents  an  eigenvalue  problem  for 
A().  However  we  defer  discussion  oT  this  problem  until  the  following 
section  (where  a  detailed  investigation  is  carried  out  of  this  aspect). 
Instead,  let  us  turn  to  consider  higher  order  terms  in  the  expansions 
(5.6)  and  (5.7).  Taking  terms 

0[h(7.)  Z  ^  c  ^  Zp  (log  Z)q'’  ) 
in  (3.4)  yields 
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MVOI)  =  A.  {-  +  WO  I  M"  -  *“']} 

♦  *  { .  man:  ♦  „,  ( ar  .sol’)}. 

(5.16) 

However,  on  account  of  (4.3)  this  equation  may  be  written  as 

L(V()|1  =  (A(  +  K()l A()>  {  ‘ 

+  Voo  (  M".  VQQ*)}.  (5.17) 

1  r  r2  11 

The  boundary  conditions  for  this  system  arc  essentially  the  same  as 
those  for  voq;  these  can  only  be  satisfied  if 

Aj  =  -  K0i  A0,  (5.18) 

implying  that 

VO 1 ( r )  =  Aoi  ¥00(r>.  (5.19) 

where  Aqj  is  a  constant  (amplitude).  It  is  straightforward  to 
determine  higher  order  terms  in  the  A  expansion,  in  a  similar  fashion. 
For  example 

a2  =  A()  K012  -  Aq  Kq2  -  i  Kqj  A0  (5.20) 

and  hence 

V02  =  a()2  V00(f).  (5.21) 

Indeed,  the  following  general  result  is  applicable 

V()n  =  A0n  VooCO-  (5.22) 

To  progress  further,  in  particular  to  determine  terms  that  arc 
0(Z'*)  smaller  than  those  considered  already,  let  us  investigate  terms 

.  M  AnZ  p-i 

0[h(Z)  Z  2  c  2  Z  (log  Z)q'‘  ] 

in  the  governing  equation.  This  yields  the  following  equation  for  yjq 
MVIO)  =  P  R1  *  Aq  R2,  (5.23) 

where 
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R,  =  ML^QCL  +  y0()  Ml  .  VQOM,,( 

r  r 

r2  =  yoo'^io'  +  Mjlrnl .  Mini’  .  (5.24) 

r  r2  r 

In  view  of  (4.4) 

t,0  =  K,„  T{)0,  (5.25) 

and  so 

R2  =  K|()R|.  (5.26) 

Repeating  the  arguments  used  to  determine  Aj  and  A 2  previously, 

then 

P  =  A()  Kl() 

=  Ja{).  (5.27) 

Finally  for  this  section,  let  us  consider  briefly  the  outer  solution, 
applicable  to  the  i)  =  0(1)  scale.  In  view  of  the  r  =  0(1)  solution, 
in  particular  its  0(log  r)  behaviour  as  r  — ♦  ■»,  together  with  (5.6), 
then  for  t]  =  0(1)  the  solution  is  expected  to  develop  in  the  following 
form 

\j/(T),Z)  =  g(Z)  h(Z)  Zp  (log  Z)(l  |vo(T),e)  +  |  Vi (rj ,e) 

+  0(Z-2)},  (5.28) 

where 


g(Z)  =  f(Z)/e.  (5.29) 

It  is  then  possible  to  obtain  an  exact  solution  for  yq  which 
matches  on  to  the  r  =  0(1)  solution.  This  is  given  by 


(5.30) 


where  Aq  is  a  constant,  and 


T0  =  I  c"  T0n(ti).  (5.31) 

n=() 


I f  we  now  expand 


(5.32) 


¥0  =  X  en  Von(n). 
n=0 

then  it  is  straightforward  to  show  that 

V00  =  A00>  (5.33) 

(which  matches  on  to  (5.15)),  and 

Voi  =  A()()  +  A01  •  (5.34) 

where  Aqj  is  an  arhitary  constant.  Other  terms  may  be  obtained  similarly. 

In  the  following  section  we  go  on  to  consider  numerical  solutions 
to  (5.10).  The  value  of  q  is  determined  in  the  Appendix. 
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The  problem  was  tackled  using  three  separate  numerical  techniques. 
The  first  comprised  a  fourth  order  Rungc-Kutta  technique,  shooting 
inwards  from  r  =  r„  (chosen  to  be  suitably  large).  In  particular, 
the  technique  involved  (i)  imposing  a  solution  of  the  form  (5.12) 
at  r<„,  generating  values  of  voqa  ( r  =  1 )  and  yoOA'  (r=l)  and  then 
(ii)  imposing  a  solution  of  the  form  (5.13)  at  rw  (or  5.14)  depending 
on  the  sign  of  Re  { i  A()i  > ) ,  generating  values  of  yooB  (r=l)  and 
V()0B'  (r=l)  (or  V()()C  (r=1)  and  Voo^'  (r=l)).  The  value  of 
Ay  was  then  chosen  by  Newton  iteration  by  imposing  (5.11),  by  forcing 
the  determinant 


vooA 

( r=l ) 

vooB  («•=*) 

vooA' 

( r=1 ) 

voo8’  (r=]> 

V00A 

(r=1) 

vooc  (r=1> 

vooA' 

(r=l) 

vooc'  («■=') 

(6.1) 


(6.2) 


to  zero. 

The  second  numerical  scheme  employed  involved  using  a  second -order 
finite-difference  approximation  to  (5.10),  constructing  a  quadra-diagonal 
system  (corresponding  to  the  approximation  to  (5.10),  together  with  the 
boundary  conditions  (5.11),  (5.15))  and  the  determinant  of  this  system  was 
forced  to  zero  by  adjusting  Aq  by  Newton  iteration. 

The  third  numerical  scheme  used  was  a  direct  (global)  finite-difference 
approach;  using  the  same  finite-difference  scheme  as  our  second  scheme, 
the  system  was  instead  written  in  the  form 


A  -  A()  B  =  0,  (6.3) 

where  A  and  B  arc  both  square  matrices.  The  Aq  were  determined 
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by  using  NAG  routine  F02GJF,  suitable  for  solving  generalised  eigenvalue 
problems  of  this  kind.  This  scheme  has  two  distinct  advantages 
(i)  of  not  requiring  iteration  and  (ii)  generating  multiple  values  (if 
present)  of  Aq  simultaneously,  however  it  can  require  substantial 
computer  storage. 

Results  from  all  three  schemes  were  found  to  agree;  in  practice  the 
procedure  was  usually  to  obtain  estimates  to  the  values  of  Aq  using  the 
third  scheme.  If  these  were  then  deemed  of  to  be  insufficient  accuracy, 
enhanced  solutions  (obtained  on  a  liner  and/or  more  extensive  grid)  were 
obtained  using  the  second  scheme  (i.c.  the  local  finite-difference 
scheme) . 

Results  were  obtained  for  a  range  of  p.  It  was  found  that  at  all 
the  values  of  p  investigated,  there  arc  many  (probably  an  infinite 
number)  values  of  Aq.  Further  all  three  methods  did  yield  a  large  number 
of  spurious  modes.  However  these  were  usually  readily  identifiable,  being 
strongly  dependent  upon  grid  size  and  range,  whilst  genuine  modes  were 
comparatively  grid  insensitive. 

Results  for  Rc{Aq}  arc  shown  in  Fig. 6  and  for  lm{A()}  in  Fig. 7. 

Just  the  first  four  modes  arc  shown  in  each  ease  -  higher  modes  become 
extremely  difficult  to  compute  (and,  indeed  distinguish  from  each  other 
and  also  the  previously  described  spurious  modes),  particularly 
in  the  limits  of  p  — »  »  and  P  — *  0.  However  the  trends  arc  clear,  namely 
that  | Af) |  — ♦  »  as  p  — *  0  and  lApi  — >  0  as  p  — »  »,  for  all  modes. 

In  the  following  section  we  investigate  these  two  limits  asymptotically. 
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In  this  section  the  limits  p  — *  »  and  p  — ♦  0  in 
equation  (5.10)  arc  considered,  for  which  some  analytic  progress  is 
possible. 

7.1.  Thfc  limit  B  — ♦  » 

Physically,  this  corresponds  to  a  low  frequency  limit  to  the 
problem.  The  numerical  results  presented  in  the  previous  section 
indicate  that  (all  the)  A()  — ♦  0  as  P  — ♦  ».  Consequently  if  r  =  0  (1), 
then  to  leading  order  (assuming  Aq  =  o(l)) 

.  m"  +  m'  =  o, 

r  r2 

with  Yoo  (1)  =  YOO'  0)  =  (7.1) 

The  solution  to  this  system  is  then 

V00  =  B0  *°K  r  -  2  +  2  } •  (7-2) 

where  Bq  is  some  arbitary  constant.  This  solution  must  ultimately 
cease  to  be  a  valid  approximation  to  (5.10)  as  r  — *  »,  specifically  when 
r  =  O(Ao'l).  Considering  the  particular  development  of  Aq  as  P  — ♦  «; 
this  is  found  to  take  on  the  following  form,  in  order  to  obtain  a  consistent 
and  meaningful  asymptotic  solution 

A0  =  Y(P)  1*0  +  e  \\  +  0(e2)|,  (7.3) 

where  e  =  -  — (7.4) 
log  y 

and  y(P)  must  be  determined  from 

y  log  (y-i)  =  p-1.  (7.5) 

This  is  a  transcendental  equation  for  the  small  parameter  y  (scc 
Duck  1984,  Duck  and  Hall  1989  for  similar  examples).  In  order  to  obtain 
a  meaningful  balance  of  terms  when  r  =  0(A(j'i),  it  is  necessary  that 
A()  =  i  (7.6) 


(the  leading  term  in  the  expansion  for  Aq). 


In  view  of  these  comments,  and  the  above  comments  regarding  the  scale  of 
r  for  which  (7.2)  ceases  to  be  a  valid  approximation  to  (5.10),  we 
define  the  outer  lengthscalc 

p  =  Y1  r  =  0(1),  (7.7) 

where  the  following  problem  must  be  considered 
Xppp  -  £  Xpp  +  Xp  [  ^2  +  i  P  +  *1  ] 

-  lx  =  {),  (7.8) 

wi  th 

X  ~  cj  |  log  p  -  i  A.j  |  as  p  -4  «, 

X  =  0  (P2)  as  p  -»  0.  (7.9) 

where  cj  is  an  arbitary  constant,  and  x  is  related  to  yoo  by 

x  =  _I__  V{,0.  (7.10) 

log  y-i 

The  system  (7.7)  -  (7.8)  represents  a  well -posed  eigenvalue  problem  for 
the  A-i's  which  was  solved  using  the  three  numerical  techniques  described 
in  the  previous  section  (indeed  (7.8)  is  very  similar  to  (5.10),  and  is 
of  about  the  same  computational  complexity,  save  for  the  absence  of  any 
physical  parameters). 

Values  for  the  first  few  Xj's  arc  tabulated  in  Table  1  (accuracy 
to  at  least  the  number  of  digits  shown).  It  appears  that  all  the 
Xj's  possessed  the  same  real  value  (and  hence  decay  rate)  to  within  the 
accuracy  of  the  computation.  The  evidence  was  that  a  large  (probably 
infinite)  number  of  these  modes  exist;  these  higher  modes  were  difficult 
to  compute  accurately,  requiring  small  grid  sizes  and  extensive  grid 
domains.  Further,  with  increasing  order,  the  imaginary  part  of  the  X] ' s 
became  progressively  more  negative,  although  the  difference  between  modes 
did  diminish.  Indeed,  these  trends  can  be  confirmed,  asymptotically, 
by  carrying  out  a  |X|i  »  1  analysis  on  (7.8)-(7.9).  In  this 
limit,  a  WKB  solution  to  (7.8)  exists  of  the  form 
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X  = 


Bjpi 


|  i | i logp+Aj |i  dp 

PO 


i logp+Xj |3/4 

+  cj| logp+iXi 1  (7.11) 

Tor  Rc  (logp-iX|)  >  0  (and  the  path  of  integration  lies  within  this 
region)  where 

P()  =  ca».  (7.12) 

(and  we  expect  ej  =  o  (B|)),  wliilst 
X  =  c | | logp-iXj | 


pi 


[ i logp+Aj | 3/4 


{*• 


|  i | i logp+Aj |i  dp 


,P0 


P 

-J  i [ i 1 ogp+X] ] i  dp 


PO 

+  A2  c 


} 


(7.13) 


for  Rc  (logp-iXj)  <  0  (and  the  path  of  integration  lies  within  this 
region) . 

A  routine  treatment  of  the  transition  layer  about  p  =  pp 
reveals 


A2  =  iAj .  (7.14) 

To  proceed  further,  consider  an  inner  layer  wherein 

p,  =  Xji  p  =  0(1),  (7.15) 

with  x  satisfying  the  following  equation  to  leading  order 

XPlPlPl  '  P  *PlPl  +  XPl(pJ  +  ~  (7-16) 

the  solution  of  which  is 

Xpi  =  j  Bppj  J0(pi).  (7.17) 

(the  second  solution  of  this  equation  involving  Yq(pi)  is  neglected 
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on  account  of  (7.9)).  Taking  the  limit  of  (7.17)  as  pj 


gives 


cos  (p,-J) 


(7.18) 


The  limit  of  (7.13)  as  p  — ♦  0  is 
X  -*  cj  | logp- iXj I 

-  }  o.») 

where 

0 

1 1  =  J  i | i logp+X| |i  dp,  .  (7.20) 

PO 

with  the  integration  path  lying  within  Re  (logp-iXj)  <  0.  If  (7.19) 
is  to  match  with  (7.18)  then 

e2I=  -1  (7.21) 

which  leads  to 

X,  =  J  -  ilog  [2  Jim],  (7.22) 

where  n  is  a  (large)  positive  integer.  For  consistency,  we  also  require 

B0  =  2i  ^  A0.  (7.23) 

Xj* 

The  formula  represented  by  (7.22)  was  used  to  obtain  asymptotic 
estimates  to  the  results  shown  in  Table  1.  Mode  II  corresponds  to 
n  =  1,  mode  III  corresponds  to  n  =  2  and  so  on;  it  is  seen  the 
agreement  between  the  compulcd  asymptotic  results  is  most  satisfactory, 

(^  =  0.785...).  It  is  quite  clear  that  this  asymptotic  form  will  fail  when 
n  =  Ofy*1). 

The  leading  order  terms,  namely  Rc(Ap)  =  ^  y  e  and  Im(A())  = 
y  arc  shown  on  Figs. 6  and  7  respectively,  for  comparison  with  the 
numerical  solutions  obtained  from  the  full  equation,  (5.10).  The 
results  arc  not  contradictary,  given  the  "largeness"  of  the  small 
parameter  e.  Indeed,  computations  for  A()  from  (5.10)  at  larger 
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values  of  p  did  become  exceedingly  difficult,  due  to  the  large 
lcngthscalc  (0(Y‘*)),  together  with  mode  "jumping"  caused  by  the 
close  proximaty  of  modes,  which  made  the  use  of  grid  refinement  with 
the  local  method  impractical. 

7.2.  The  limit  B  —>  0 

This  corresponds  to  the  high  frequency  limit  of  the  problem. 

According  to  the  numerical  results  presented  in  Section  6  tAqi  increases 
as  p  — »  0.  This  limit  is  now  investigated. 

It  is  possible  to  write  a  WKB-typc  approximate  solution  to  (5.10) 
(assuming  iAqi  and  p*1  arc  both  large)  as 
Y00<r)  =  Bl  [l°8  r  -  ] 

r  i  r5/4 

+  |A{)  log  r  -  p  J  ri 


B2 


Iaq  tog  r  -  rr  jj  ar 


0 


-iffA()  log  r  -  ill  dr  i 
+  B3  e  1 


(7.24) 


(for  Re  {A()logr  -  ^  }  <  0,  and  the  path  of  integration  lies  within 
this  region),  where 


.  -m 

B3 - |  Y  - 

1  1  '  P  )  +  A° 

-1/4 


3/2 


B1  =  *  A()  (-  P  )  [  B2  c'  +  b3  c''  ]  - 

and  ro  -  «p  [  ] 

is  the  turning  point,  and 

i  ii 


r() 


I  =  -  f  [A0  l08  r  -  p  ]  dr. 

1 


(7.25) 

(7.26) 

(7.27) 


(7.28) 
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(7.25)  and  (7.26)  arc  obtained  by  imposing  boundary  conditions  on  r  =  1 , 
and  the  integration  path  lies  within  Re  {A()logr  -  p)  <  0. 

For  Re  (A()logr-  >  0  the  WKB-lype  approximate  solution  can  be 


wri t  ten 


¥()<>('>  =  »l  [lo8  r  *  ] 


{()[  i  ■  A<>  lo» r  ] 


[  jj  -  A()  log  r  ] 


.  (7.29) 


where  it  has  been  assumed  Rc{(-Aq)J}  <  0  (otherwise  we  require  the 
negative  root  inside  the  integral),  B|  is  given  by  (7.26),  and  the 
integation  path  lies  within  Re  {Aqlogr  ‘  p)  >  0.  In  order  that 
(7.29)  matches  to  (7.24)  across  the  transition  layer  of  thickness 
0(Aq’  ),  (routine)  treatment  (see  also  the  analysis  for  Aq  =  0(P  '  ) 


below)  of  the  latter  yields 
B|  =  -  i B2 , 


(7.30) 


and  the  fol lowing  di spers ion  relationship  for  Aq  results 


:  J/2 


:  ,3/2 


(  -  P  )  -  i  Ao  =  c  {i  (  -  p  )  -  A(,  }.  (7.31) 

It  turns  out  that  there  arc  two  distinct  families  of  solution  as 


p  -4  0. 


The  first  family  of  solutions  as  P  — »  0  corresponds  to  Aq  =  0(P"') 


More  spcci fical ly 


A()  =  P'1  [  A0  +  P1/3  A|  +■  •  •  ]  • 


(7.32) 


where  Aq ,  Ai  arc  generally  0(1)  quantities.  This  implies  that 
tq-1  =  0(1).  Consequently  to  leading  order  (7.31)  reduces  to 


2i  I 

c  =  - 1 . 


(7.33) 


However  it  appears  1  =  0(P'^),  and  so  there  is  a  contradiction, 


which  can  only  be  avoided  if 
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(7.34) 


r0 

|  |Aq  log  r  -  i 1*  dr  =  0, 

1 

or 

i 

A() 

|  c*p  p*  dp  =  0,  (7.35) 

0 

or 

v  (  T  f  1  -  (7  M> 

A<) 

where  y(zl*z2)  represents  the  incomplete  gamma  function.  This  represents 

A 

an  eigenvalue  problem  for  Ap,  which  was  solved  numerically  using  a 
combination  of  trapezoidal  quadrature  and  Newton  iteration;  results 

A 

for  the  first  few  Ap  arc  shown  in  Table  2.  Note  that  there  appear 
to  be  many  values  (probably  an  infinite  number),  although  these 
seem  to  be  concentrated  within  a  finite  annular  region  in  the  complex 

A 

A()  plane.  As  the  order  increased,  the  values  become. very  close  to 
neighbouring  values,  and  the  computation  became  exceedingly  difficult; 

A 

however,  with  increasing  order  the  values  of  Ap  do  seem  to  be 
approaching  a  finite  value  (indeed  the  author  was  unable  to  find 

A 

any  solution  for  tApi  <  0.098).  Note  too  that  it  is  easy  to  show 
using  integration  by  parts  that  there  arc  no  solutions  to  (7.35)  as 

A 

|A()|  — »  «,  whilst  using  the  asymptotic  expansion  for  the  incomplete 
gamma  function  (Abramowitz  and  Stegun  1964)  it  is  also  possible  to  show 
that  no  solutions  exist  as  i Aq i  — »  0  either;  this  then  confirms  our 
statement  about  the  values  of  Aq  being  confined  to  an  annular  region  in 
complex  Aq  space. 

A  A 

Note  also  that  both  Aq  and  -complex  conjugate  {Ap}  are  roots  of 
(7.35);  however  the  latter  family  of  solutions  may  be  disregarded  since 
in  all  eases  we  require  Aq  to  possess  a  positive  real  part. 
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The  second  family  of  solutions  for  Aq  occurs  when 
Aq  =  0(p'3/2).  in  this  ease,  from  (7.27),  trQ-li  «  1, 
and  indeed  the  wall  (r=l)  lies  inside  the  transition  layer. 
Consequently,  we  arc  unable  to  use  (7.31),  but  must  consider  the 
transition  layer  in  detail  (although  this  is  quite  a  routine  task). 
Suppose 

A(,  =  P‘3/2  A0,  (7.37) 

where  Aq  =  0(1).  Then  defining 

C  =  ( r- 1 )  P'1.  (7.38) 

to  leading  order  (5.M))  reduces,  to 

AAA  A  A 

vooccc  +  Vooc  (AoC-i) 

-  A0  Voo  =0.  (7.39) 

Writing  1  =  (-A0)'1/3  o  +  i-  ,  (7.40) 

Ao 

A 

and  differentiating  (7.39)  with  respect  to  C,  yields 


^^oooo  '  0  ~ 

The  required  solution  (that  is  not  exponentially  large  as 
i  s 


(7.41) 
o  — »  «°) 


Y()00  =  0  Ai  (o) ,  (7.42) 

(where  D  is  independent  of  a). 

The  implementation  of  the  boundary  conditions  on  £  =  0  requires 
Voo£££  (C=0)  =  o,  and  so 

~  1  /3 

Ai'  (  -ULM -  ]  =  o.  (7.43) 

A() 

Now  since  the  zeroes  of  the  Airy  function  and  its  derivative  arc 
confined  exclusively  to  the  negative  real  axis,  then 

Ai*  (-Cn>  =  0.  (n=l ,2,3, ... )  (7.44) 

where  the  £n  arc  real  and  positive  and  tabulated  by  Abramowitz  and 
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Stcgun  (1964).  Consequently 

s0  -  -77-3/2+  W''  <7«> 

where  the  appropriate  roots  have  been  chosen  to  ensure  boundedness 
of  the  Airy  function. 

It  is  interesting  (although,  in  some  ways  not  too  surprising)  that 
(7.45)  is  identical  to  the  corresponding  expression  found  in  the 
analogous  planar  study  (Lam  and  Roll  1960,  Ackcrbcrg  and  Phillips 
1972  and  Goldstein  1983),  although  of  course  the  corresponding  f (Z) 
is  quite  different  in  the  present  ease. 

As  a  check  on  the  numerical  results  as  p  — ♦  0,  on  Pig  8  the 
variation  of  p 3/2  will,  p  js  shown  (first  three  modes).  It  is  very 

clear  that  these  results  approach  those  given  by  (7.45)  as  p  — »  0. 

The  0(P_1)  family  of  results  refer  to  higher  modes,  and  thus  it  is 
not  realistically  possible  to  compare  our  numerical  results  with  this 
fami ly . 

In  the  following  section  we  draw  some  conclusions  from  this  work. 


1L  Conclusion 


In  this  paper  the  cfTcct  of  smalt  amplitude  frccstrcam  oscillations 
on  an  otherwise  steady  boundary  layer  on  an  axisymmetric  body  has  been 
investigated.  Particular  attention  has  been  focused  on  the  far- 
downstream  eigenvalues  and  cigcnsolut ions.  As  noted  in  Section  1,  in 
the  ease  of  the  planar  problem,  two  distinct  families  of  cigcnsolut ions  have 
been  presented,  namely  those  originally  considered  by  Lam  and  Rott  (I960) 
and  those  considered  by  Brown  and  Stewartson  (1973a, b),  with  the  former 
family  having  decay  rates  that  decrease  with  increasing  order,  whilst  the 
latter  family  have  decay  rates  that  increase  with  increasing  order. 

In  the  present  study,  eigenvalues  appear  to  occur  with  decreasing  decay 
rale  with  increasing  order.  However,  some  of  the  asymptotic  work  in 
Section  7  (in  particular  that  relevant  to  P  — ♦  »,  with  Aq  =  0(P‘l)) 
docs  strongly  suggest  that  a  finite  value  of  Aq  is  being  approached 
with  increasing  order.  Indeed,  the  author  was  unable  to  obtain  a 
consistent  asymptotic  solution  to  (5.10)  for  p  =  0(1),  Aq  — ♦  0,  again 
suggesting  the  finite  limit  of  Aq  with  increasing  order.  This,  in 
some  ways  may  be  regarded  as  a  rather  more  satisfactory  state  of 
affairs  than  that  found  with  the  Lam  and  Roll  (1960)  cigcnsolut ions , 
which  have  decay  rates  that  become  dimini  singly  small  with  increasing  order 
(although  see  our  comments,  attributed  to  Goldstein  ct  al  1983,  in 
Section  1).  Further  the  P  — *  0  work  of  Section  7  docs  suggest 
that  all  modes  possess  the  same  decay  rate  in  this  limit  up  to  at  least 
second  order. 

However,  it  may  well  be  that  the  planar  work  of  Brown  and  Stewartson 
(1973a, b)  could  perhaps  be  extended  to  include  the  effects  of  curvature,  to 
yield  a  further  (perhaps  related)  family  of  cigcnsolut ions .  A  further 
interesting  study  would  be  an  investigation  of  the  far-downstream  evolution 
of  the  cigcnsolut ions.  Just  as  in  the  planar  ease,  these  all  become 
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increasingly  oscillatory  far  downstream,  and  will,  as  a  consequence, 
ultimately  cease  to  be  valid  approximations  to  the  Navicr  Stokes  equations. 
This  will  lead,  presumably,  to  the  formation  of  unstable  Tollmicn- 
Schlichting  waves,  in  a  manner  analogous  to  that  described  by  Goldstein 
(1983)  in  the  planar  ease. 

However,  there  arc  a  number  of  (other)  important  differences  between  the 
planar  and  the  axi symmetric  cigcnsolut ions  and  eigenvalues.  Most 
importantly  the  downstream  (i.c.  axial)  behaviour  of  these  cigcnsolutions 
(described  by  l'(Z))  which  is  quite  different  in  the  two  eases,  in  the 
axisymmctric  ease  being  given  by  (5.6)  whilst  in  the  planar  Blasius  ease 


(8.1) 


as  shown  by  Goldstein  (1983),  (where  x  is  the  slrcamwisc  coordinate). 

Note  that  if  the  basic  flow  were  of  the  form  T  =  xm  F(q),  with  q  =  y/xm,  y 
being  the  transverse  boundary  layer  variable,  then  using  arguments  similar 
to  this  paper, 

f(Z)  =  xp  cxp{  -A  x2ll-m+l},  for  2n-.n+l  >  0.  (8.2) 
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Appendix 


Consider  terms 

*  M  M 

o{h(Z)  Z  *  c  *  Z  P_1  (logZ)  q'1  J 

in  (5.1)  to  enable  as  to  determine  the  value  of  q. 

It  is  fouiut,  after  some  algebra 

l{¥1!}  =  P{  mim’  +  yoi'm' 

■  *0  [  *1"  •  M  ]  .*,[  SB"  .  SB’]} 

r  r2  r  r2 

+  q{M3a'  .  «  ^  (Stan",  an'  ]} 

+  a,  {  fon'vio'  .  Imlsm'+  yo  [  Im"  .  !u’  ) 

1  r  r  1  r  r2  ' 

+  ,lo[y:M']} 

+  a0  {-  .  hii^rni .  hi'm: 

+  xf(){)  j  In  .  Ill  j  +vo,  [  Iiii  -  Ilq  ] 

r  r2  r  r2 

r  r2 

*  —  ¥<)()'  %)’  +  j  %)  [  ^  -  ^1  • 

2r  1  1  r  r2j 

(A.  1 ) 

Recalling  the  expressions  obtained  already  for  p,  Aj ,  vqj ,  vjq. 

(A.l)  leads  to  the  following  (slightly  simplified)  equation 

Mfll)  *  0(r)  {  A()  K]()  Aol  +  ■}  +  A()  K()|  A[0  -  <A|(>  +  K , (,)  A„ 

.  JlliBl  ♦  ¥0fl  [  ill  '  .  Ill'  ] 

r  1  r  r2  J 


where 
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(A. 3) 


Q(r)  =  WQQ1  ]. 

r  1  r  ■* 

It  is  quite  clear  that  \|tj|  =  0(r2(log  r)2)  as  r  — ♦  To 

simplify  arguments  later,  we  write 

Yll  =  V1 1 1 ( r)  +  q  VI 1 1 1 < r) ,  (A. 4) 


where  \|fllH(r)  is  any  regular  function  which  has  the  following  behaviour 
as  r— »  oo 

Vll,,(r)  =  -  p)^2  +  o(D.  (A. 5) 

The  governing  equation  for  \|/jjl  may  then  be  written  in  the  form 
L(Vl  1 1 1  =  R3<r)  +  <1  r4  («■)*  (A. 6) 


where 


R3(r)  =  Q(r)  {  A()  K10  A0j  +  A0  k()]  A10 

*  (Alt)  +  R1())A()}  -  -1 1-^00- 


and 

R4(r)  =  -  j  Q(r)  -  L{Vl,"}.  (A. 7) 

The  boundary  conditions  that  must  be  applied  to  this  system  arc 

Vll'O)  =  -  V|  1 1 1  ( 1 ) 

Vl  1 1  *  (1 )  =  -  Vll11 ' 0) 

and  that  yj|l  — ♦  0  as  r  — »  oo,  (A. 8) 

The  value  of  q  is  then  determined  by  the  condition  that  a  solution 
to  this  system  exists. 

Consider  now  the  (complex  conjugate  of  the)  adjoint  to  the 
system  (5.10),  denoted  by  ^(r),  and  determined  from 
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(A. 9) 


V+”  '  +  t*”  +  V+’  [  A0  log  r  -  i  -  1  1 
1  P  r2J 

+  2  Jlf+rO, 

subject  to  the  boundary  conditions 

V+(r=l)  =0,  y+,(r=l)  =  1  (say),  (A.  10) 

=  o  ((logr)*2)  as  r  — »  «.  (A. 11) 

If  we  now  further  suppose,  as  we  arc  quite  at  libery  to  so  do  (although 
this  simplifies,  but  is  not  crucial  for  our  arguments)  that 

yll  H(i)  =  yll  n*(l)  =  0,  (A.  12) 

then 

oo 

|  R3(r)  y+(r)  dr 

q  =  -  1 -  (A. 13) 

«4(r)  V+(r)  dr 

(Hartman  1964,  for  example). 

This  (at  least  in  principle)  determines,  or  provides  a  means  of 
determinining  the  index  of  the  logarithmic  term  mul t iplying  f(Z). 

At  this  stage  it  would  also  appear  to  be  legitimate  to  set  the  function 
h(Z)  in  (5.6)  equal  to  a  constant,  although  categorical  determination 
of  this  point  seems  difficult  because  of  the  algebraic  complexity 
in  extending  the  analysis  to  higher  order. 
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Mode 

*1 

Asymptotic  Xj 

I 

.785  +  . I60i 

1! 

.785  -  1 . 282 i 

J  -  1 . 266 i 

Ill 

.785  -  1.9341 

J  -  1.9591 

IV 

.785  -  2 . 340 i 

■1 

.785  -  2.631 i 

J  -  2.652i 

VI 

.785  -  2.858' 

J  -  2.8751 

VII 

.785  -  3 . 044 i 

J  -  3 . 05  7  i 

VIII 

.785  -  3.2()()i 

J  -  3.21 1 i 

Tabic  1  Values  of  Xj 


A0 

.1408 

+  .2262X10' * i 

.7455x10-* 

+  . 7991X10-2 i 

.5076XI0'1 

+  .4181X10-2j 

.3848X10** 

+  .2603X10" * i 

.3099X10-* 

+  .1790X10-2) 

.2594X10-* 

+  .1313X10-2i 

.2231X10** 

+  .1007X10-2) 

.1952X10-* 

+  .7996X10'3! 

.1742X10'* 

+  . 65 14X10-3  j 

.157IXIO-* 

+  .54I8X10-31 

A 

Table  2  Values  of  Aq 
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^rrjr  =  1 


Fig.  1.  Variation  of  ¥orrlr-i 


Yrr|r  =  1 
A 

2.0 


y  /~  Re{\}/rr|rsi} _ Asymptotic 


1.5 

1.0 


Asymptotic 


MVrr|r  =  l} 


_J _ i _ i _ i _ i _ i _ i _ i - i - 1 _ i _ 

1  23456789  10  11 


Fig.  2.  Variation  of  Vrrir_,.  P  =  i 
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0  123456789  10  11  z 


Fig.  3.  Variation  of  vr r I r_ j>  P  =  1 


Vrr|r  =  1 


0  123456789  10  11z 


Fig.  4.  Variation  of  yrr|  j,  P  =  2 
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